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We show theoretically that a monopole defect, analogous to the Dirac magnetic monopole, may 
exist as the ground state of a dilute spin-1 Bose-Einstein condensate. The ground-state monopole 
is not attached to a single semi-infinite Dirac string, but forms a point where the circulation of 
a single vortex line is reversed. Furthermore, the three-dimensional dynamics of this monopole 
defect are studied after the magnetic field pinning the monopole is removed and the emergence of 
antimonopoles is observed. Our scheme is experimentally realizable with the present-day state of 
the art. 



Introduction. — Topologically nontrivial configurations 
of quantum and classical fields play a fundamental role 
in the physics of various phase transitions ranging from 
superfluids to the early universe [ll,[2[. In particular, they 
can give rise to exotic quantum states relevant to elec- 
tromagnetism (sj, elementaryparticles grand unified 
theories [H , and cosmology ]6| . Given the difficulties of 
making detailed observations on cosmological scales or 
probing conventional matter on ultrashort length scales, 
experimental evidence for many of these topological ex- 
citations is scarce and indirect at most. 

Dilute ultracold Bose gases with a spin degree of 
freedom can host a variety of topologically interesting 
structures, such as coreless vortices, knotted textures, 
skyrmions, and several types of monopoles The 
high controllability of these quantum systems and their 
observability using optical imaging techniques allow, in 
principle, detailed experimental investigations of the vari- 
ous defects. In practice, however, topologically nontrivial 
configurations occur as excited states and their topolog- 
ical nature dictates that they cannot be easily created 
from a topologically trivial ground state. In particular, 
configurations relevant to tabletop experiments of exotic 
phenomena in high-energy physics and cosmology can be 
especially challenging to create in an atomic gas since 
they require intricate manipulation of the atomic cloud 
in all three spatial dimensions [l6l-[2Qj. 

A major step forward was taken recently when a robust 
method for creating an analogue of the Dirac monopole 
was proposed [2o[. The Dirac monopole [3] is the sim- 
plest model for a magnetic point charge which unlike its 
electric counterpart, has not been convincingly observed. 
In Ref. [20I, the monopole defect was created into the 
spin texture of a dilute Bose-Einstein condensate (BEC) 
by adiabatically modifying external magnetic fields. The 
resulting monopole state is not the ground state of the 
system because of its tendency to degrade due to dynam- 
ical instabilities associated with the Dirac string. In this 
Letter, we demonstrate that an analogue of the Dirac 
monopole in an atomic gas can exist as a ground state 
configuration. We find that in an experimentally feasible 
time-independent magnetic field the ground state of the 



BEC corresponds to a strong-field seeking state (SPSS) 
with a monopole defect. A similar defect is also found 
for the weak- field seeking state (WFSS). In contrast to 
the previous studies, we find that the minimum-energy 
monopole is not attached to a single Dirac string with two 
circulation quanta but manifests itself as a point where 
the circulation of a single-quantum vortex is reversed. 
A similar configuration may be obtained from the ideal 
Dirac monopole by applying a gauge transformation [2l| . 
Our results show that the Dirac monopole is a physically 
viable concept and can represent a robust and long-lived 
state. 

We consider a condensate of ^^Rb atoms with ferro- 
magnetic spin-spin interactions [22|. The total hyper- 
fine spin of the constituent atoms is F = 1, and the 
order parameter is a three-component spinor field. In 
the presence of a strong external magnetic field, the con- 
densate spin tends to align with the local field. Hence, 
external fields can be used to imprint topologically non- 
trivial spin textures to the condensate. For topological 
reasons, the ferromagnetic phase cannot sustain stable 
isolated point-like defects [23|. This does not, however, 
completely exclude the existence of point defects, as they 
may appear as endpoints of quantized vortices. These 
defects are characterised by a geometric charge Q^d cor- 
responding to the area on the unit sphere covered by the 
condensate spin as one encloses the point defect in the 
spatial coordinate space (2q| . 

Let us consider an external magnetic field which is a 
combination of two crossing quadrupole fields. 



(1) 



Since S is a monopole-free field. Maxwell's equation 
V • -B = imposes a condition 2B[ -\- B2 = 0. The 
alignment of the hyperfine spin with the magnetic field 
in Eq. ([T]) gives rise to a spin texture that is depicted in 
Fig. [TJa) and has the geometric charge Qsd = 1- The 
analogy to the Dirac monopole comes from the vortic- 
ity fig = V X i;^, where Vs is the superfluid velocity. 
The vorticity fls is equivalent to the magnetic field of a 
magnetic monopole and the Dirac monopole can be con- 
sidered as a point source of the superfluid flow. Using 
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the Mermin-Ho relation [2^, the geometric charge Q2,d 
of the spin texture can be related to the total vorticity in 
the system. Vorticity corresponding to the spin texture 
in Fig.flja) takes almost everywhere the radially outward 
hedgehog form 
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in the scaled units {x' ^ y\ z') = (x, 2z) This estab- 
lishes the analogy to the magnetic monopole proposed by 
Dirac 

Mean- field theory. — In the zero-temperature limit, the 
stationary states of the condensate are solutions tojthe 
time-independent Gross-Pitaevskii (GP) equation (25| 

= /i^(r), (3) 

where the Hamiltonian for a spin-1 condensate is 

= h{r) + col^ (r) p + C2^\r)J^^{r) • JT. (4) 

Here, !F = (J^^, J^y, J^z)^ is a vector of spin-1 matri- 
ces and h{r) is the single-particle Hamiltonian given by 
h{r) = -Ti?V'^ /2m^Vopt(r)^gFlJiBB{r)'T, where gp is 
the Lande ^-factor and /i^ is the Bohr magneton. Unless 
otherwise mentioned, we consider a three-dimensional 
optical trapping potential V^pt(^) = muj'^r^ /2. The cou- 
pling constants are given by cq = 47rfi^(ao + 2a2)/3m and 
C2 = 47rfi^(a2 — ao)/3m, where aj is the 5- wave scattering 
length corresponding to the channel with total hyperfine 
spin /. The Hamiltonian in Eq. (ffl corresponds to a 
free-energy functional of the form [25|, [26| 



Em = jdr 



|-|V^(r)|V[Fopt(r)-/i]|^(r)|^ 



^\S{r)\^^gF^BB{r)-S{r) 
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where S = ^"f"^^. The dynamics of the BEG are solved 
from the time-dependent GP equation, iTidt^ = 
where 1-L is given in Eq. (|4]). 

Let us write the stationary order parameter in the 
form ^ = where is a scalar field and is a unit- 
normalized spinor field. Assuming that the condensate 
spins align with the local field, the free energy of the 
system can be written as (27| 
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where V is the interaction potential. Equation (|6]) is a 
Hamiltonian for a scalar particle with order parameter 
and it is equivalent to the Hamiltonian for charged 



particles coupled to a vector potential A. Hence, the 
equation of motion for the condensate is equivalent to 
that of charged particles in an electromagnetic field. The 
vector potential of the electromagnetic field corresponds 
to the superfiuid velocity of the condensate if we set 
mVs = qA [27|. 

For an optically confined BEG, minimization of the en- 
ergy in Eq. (j5j) results in a condensate in the SPSS which 
is lower in energy than the WFSS. We stress that only 
the SPSS monopole represents a robust ground-state con- 
figuration. In order to study a realistic minimum-energy 
WPSS, we consider a purely magnetic trap provided by 
the two quadrupole fields in Eq. ([1]). 
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FIG. 1. (a) Spin texture corresponding to the magnetic field 
in Eq. ([T]). (b) Vorticity fig corresponding to the spin texture 
in panel (a). Vector fields in panels (a) and (b) are symmetric 
with respect to rotations about the axis depicted with the 
dashed line, (c) Illustration of the polar angles ot and P used 
to parametrize the local magnetic field B. 

We parametrize the magnetic field using polar 
angles a and /3 [Fig. Hl^c)] such that B (r) = 
1^1 (sin /3 cos a, sin /3 sin a, cos /3) . As the Lande g- 
factor is negative for spin-1 ^^Rb and we are consider- 
ing the WFSS, the hyperfine spin S is antiparallel with 
the local magnetic field. In the eigenbasis of the spin 
operator F^, the order parameter can be expressed as 



e-*" sin^ (;9/2) \ 

sin/3 
e^" cos2 (/3/2) / 



/e-'X, (7) 



where / is the amplitude of the order parameter and C, 
if fixed by Eq. ([T]) . We substitute the order parameter in 
Eq. ([7]) into the stationary GP equation (j3j) and denote 
^ = /0C7 where /e~*^ = /q. This gives a reduced GP 
equation which can be written in the dimensionless form 



i[v2 + 2(c^vc)-v+(c^v2c)^ 

s|/o + (co + C2)/o^ 

[iBBjTlUJr, fo 



fo 
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Here ji = fi/TiuJr^ B 
CiN / (a'^hujr) and = ^^h/ {mujr)- We normalize the 
complex- valued function fo to unity as J|/opdf = 1. 



The minimum- energy density and phase distributions 
are solved from Eq. ^ using the standard imaginary- 
time propagation combined with finite- difference meth- 
ods. The strong-field seeking ground state is found by 
minimizing the free energy in Eq. (|5j) with full spin 
degrees of freedom using the successive overrelaxation 
method. The temporal evolution of the monopole defect 
is solved from the time-dependent GP equation using the 
split-operator and Crank-Nicolson methods. 

Numerical results. — The mass of a ^^Rb atom is taken 
to be m = 1.44 x 10~^^ kg, the Lande g factor is gp = 
— 1/2 and the total number of atoms TV = 8 x 10^. For 
the coupling constants of ^^Rb, we use the value C2/C0 = 
-0.01 and adopt cq = 7500. For ujr = x 250 Hz, 
the dimensional values for the parameters are given by 
B[ = 0.1 T/m for the simulations with the WFSS and 
B[ = -0.03 T/m with the SFSS. The volume considered 
in the simulation is 23 x 23 x 28 in the units of af. The 
computational grid consists of 141 x 141 x 161 points. 

For both SFSS and WFSS, the minimum-energy con- 
figuration corresponds to a monopole defect associated 
with two vortices (Dirac strings) each carrying a sin- 
gle quantum of angular momentum. Both strings carry 
vorticity towards the monopole defect which lies in the 
zero point of the magnetic field. In Ref. the Dirac 
monopole is associated with only one vortex line that 
carries two quanta of vorticity. For energetic reasons, it 
is natural that this two-quantum vortex splits into two 
singl y q uantized vortices in the ground-state configura- 



the hedgehog form of Fig.fTJb) near the monopole defect. 
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In the WFSS, the Dirac strings lie on the z-axis 
because the magnetic trap is strongest along this di- 
rection. For the SFSS, the confinement in the z- 
direction is weaker than in the x- and ^-directions due 
to the symmetric optical trap, and the two strings 
lie in the x?/-plane. For the SFSS, we also consider 
an asymmetric optical trap of the form V^pt (r) = 
m {cjO^x'^ + ^yV'^ + ^?^^) general, the Dirac strings 

associated with the monopole defect lie along the direc- 
tion that minimizes the length of the vortex lines. Parti- 
cle densities for the SFSS and WFSS are shown in Fig.[2l 
For the WFSS, the vortices corresponding to the Dirac 
strings are singular and the particle density vanishes at 
the vortex core. In the case of the SFSS, particle density 
is only partially depleted along the two vortices, implying 
that they are polar-core vortices [30] . The particle densi- 
ties for each spin state in the WFSS and SFSS are shown 
in Fig. m In the SFSS, the two vortices are manifested 
as depletion in the spin density, see Fig. ID^a). 

Vorticities and spin textures corresponding to SFSS 
and WFSS were found to be qualitatively similar, and 
we only present further results for the SFSS. The con- 
densate spin is shown in Fig. ID^b) and one observes that 
the spin indeed tends to align with the local magnetic 
field. From Figs. [5fa)-(c) we observe that vorticity takes 
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FIG. 2. Particle densities corresponding to (a) WFSS and 
(b) SFSS for X > 0. The densities are given in the units of 
N/a^. The field of view is (a) 8 x 16 x 8 and (b) 6 x 12 x 14 
in the units of a^. 




FIG. 3. Densities of the three spinor components for (a) 
WFSS and (b) SFSS. In both panels, the top segment cor- 
responds to 1^1 1^, the central segment to |^o|^, and the bot- 
tom segment to The colorbar scales are in the units 
of N/a^. The segments are bounded by density isosurfaces 
for spinor components for x > 0. On the plane x = 0, the 
isosurfaces are capped with a density colormap for the cor- 
responding spinor component. In panel (a), all isosurfaces 
correspond to density 1.4 x 10~"^A^/a^ and the singly quan- 
tized vortices are manifested as density depletion along the 
z-duxis. In panel (b), the vortices lie on the y-axis and the 
isosurfaces correspond to densities 1.2 x 10~^ N/a^ for |^i|^ 
and 1^-1 1^ and 4.5 x IQ-^AT/a? for |^o|^ 

Next we study the dynamics of the monopole after the 
external magnetic fields are turned off. We take the ini- 
tial state to be the strong-field seeking ground state with 
vortices on the y-ax.is. The fields are turned off instan- 
taneously and the state is evolved up to lO/o;^ with a 
timestep 10~^/cj^. During the time evolution, the Dirac 
strings are observed to expand, see Figs. [5jd)-(f). Fur- 
thermore, two antimonopoles emerge from the bound- 
ary of the condensate and move towards the center of 
the trap. The emergent antimonopoles are character- 
ized by Dirac strings that carry vorticity outwards from 
the monopole. Similarly to the ground-state monopole, 
the Dirac strings are coreless vortices. We note that 
the Dirac strings associated with antimonopoles guar- 
antee that the total vorticity remains zero after the an- 
timonopoles emerge. A schematic representation of the 
temporal evolution is shown in Fig. [6] and high resolution 
figures of the simulation are available in the supporting 
online material ^j. At the end of the simulation, the 
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FIG. 4. (a) Spin density for the strong- field seeking ground- 
state monopole in the z — plane. Spin density is depleted 
along the two vortices, which lie on the ^/-axis. The fig- 
ure corresponds to spin densities from 0.1 x 10~^A^/a^ to 
0.92 X 10~^A^/a^. (b) Spin of the ground-state monopole. 
The arrows represent the projection of the spin to a plane 
corresponding to x = 0. The Sx component is zero in this 
plane. 



\ \ \ \ \\ \\ 
\\\\\\\\\ntttit 
t f f f / / M t M \ \ 1 1 t 
r,v\\\\\ \ \ \ / // // // 
ZsNNVWN \ / f / / y ^ ^ ^ 



M 



'///// / I \ \ www 
/////// i \ \ \ \ \\ \\ 
/ / / / / / / I 1 \ \ \ \ \ \ \ 
\ \ \ \ \ \ \ 1 I / / / / / / / , 



-2 a 2 

x/ar 
(d) t = 2/uj,. 



7 r f r f ! t J 

\\\\ \ \ \ \ 


1 1// \ A/ \ ' 

If////// 


t t f f / / / f 


t t \ U M t 


\ \ \ \ 


f / ////// 


wwv\ \ \ 


f / / / /• y ^ y' 




t / y^^^^^ 




1 y^. 




\ \ 


'^^^z / / I 




' y ^ z / / / / 


I \ \ WWN 


''//Z / / / / 


I \ \ \ WW 


1 u n t i \ 


I I \ \ \ \ \ \ 


\\\\\\\\ 


\ n 1 1 1 1 1 



(c) 



x/ar 



(e)j = 6/ujr 



<.():; 
— 2^ 
-4 



t t t t f m 



" PM M 

\\\\\\\\\\iiftn 

,K\\\\\NM / // // // 
\\ \ f / / yyyy 



■^■^''^ / I \ \ 

y / / / / J \ \ \ WXN 

'//// / J I \ \ \ 

/////// j M \ \ \ \ W 

\ \ \ \ \ \ W M / / / / /J 



4 



-2 0^ 2 
x/ar 



(f) t = 8/uJr 

2i! 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 



FIG. 5. (a)-(c) Vorticity of the ground-state monopole at 
different locations, (d)-(f) Temporal evolution of vorticity at 
fixed location. In panels (a)-(c), arrows represent a projection 
of the unit vorticity Cls = fts/\^s\ to x^-plane for (a) y = 
0.8 X ttr, (b) ^ = X ttr and (c) y = —0.8 x ar- In panels (d)- 
(f) y = —0.7 X ttr at time instants (d) t = (e) t — Q/uJr 

and (f) t — S/uJr- The y component is presented with the 
color map. 



vorticity of the monopole defect has spread out and is no 
longer strictly of the hedgehog form. 

In conclusion, we have shown that a monopole defect 
may exist as the ground state of a dilute ferromagnetic 
spin-1 BEG in the presence of an experimentally feasi- 
ble magnetic field configuration. Vorticity of this de- 
fect is analogous to the magnetic field of a magnetic 
monopole. The ground-state monopole is associated with 
two Dirac strings that extend symmetrically outwards 
from the monopole defect. From the relation between 
the geometric charge of the monopole defect and the to- 





FIG. 6. Schematic representation of the temporal evolution 
of the monopole defect. The spheres represent the atomic 
cloud and the black lines correspond to vortices. The arrows 
illustrate the direction of vorticity associated with the single- 
quantum vortices. The condensate is depicted at two instants 
of time t = 2/ojr and t = 10/oJr- In the top panels the con- 
densate is viewed along the ^/-axis and in the bottom panels, 
along the z-axis. 



tal vorticity in the system, it follows that each of the 
two coreless vortices associated with the SFSS monopole 
defect carry a single quantum of angular momentum. 
This is remarkable, since in general the angular momen- 
tum associated with coreless vortices is not quantized. 
Thus, the monopole defect is responsible for the quan- 
tization of coreless vortices similarly to the Dirac mag- 
netic monopole which imposes the quantization of elec- 
tric charge. If monopoles are to appear spontaneously in 
BECs, one would expect for energetic reasons that they 
are associated with two singly quantized vortices instead 
of a single two-quantum vortex. Hence such monopoles 
are of great interest. In the experiments, the ground- 
state monopole can be created during the cooling process 
in situ, without the need to adjust the external magnetic 
fields in time. These monopoles are expected to be ex- 
tremely robust and long lived. 
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